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Abstract. We calculate the spin-drag relaxation rate for a two-component ultracold 
atomic Fermi gas with positive scattering length between the two spin components. 
In one dimension we find that it vanishes linearly with temperature. In three 
dimensions the spin-drag relaxation rate vanishes quadratically with temperature for 
sufficiently weak interactions. This quadratic temperature dependence is present, 
up to logarithmic corrections, in the two-dimensional case as well. For stronger 
interaction the system exhibits a Stoner ferromagnetic phase transition in two and 
three dimensions. We show that the spin-drag relaxation rate is enhanced by spin 
fluctuations as the temperature approaches the critical temperature of this transition 
from above. 
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1. Introduction 

Interest in electronic transport ranges from everyday applications to fundamental 
physics. One of the most interesting phenomena that spans this entire range, is the 
influence of a thermodynamic phase transition on the electrical conductivity. The most 
direct example is the phase transition from a normal conductor to a superconductor 
characterized by a vanishing resistance. The applications of this phenomenon are 
ubiquitous and the basic physics that underlies the transition in superconductors, the 
Bose-Einstein condensation of fermionic pairs, has emerged in research flelds from 
astroparticle physics [1] to cold- atom systems [2^ 3j. 

A system in between the latter two temperature extremes, in which analogies 
of superconductivity have been predicted, is that of a two-dimensional electron-hole 
bilayer [4^ [5]. In this case the pairs that condense are excitons formed by electrons 
from one layer with holes in the other. The relevant transport probe is in this case the 
Coulomb drag measurement [6j: a current / is driven through one layer, known as the 
''active" layer, causing a voltage drop in the other. As the layers are separated 
by an essentially impenetrable tunnel barrier, the voltage drop is predominantly 
caused by Coulomb scattering, and the drag resistivity Pd = ^d/^ has the, up 
to logarithmic corrections, characteristic quadratic Fermi-liquid-like low-temperature 
dependence pj^^ (xT^ . When the excitons undergo Bose-Einstein condensation, however, 
the drag resistivity is predicted to jump from the relatively small value proportional to 
to a value equal to the ordinary resistivity of the active layer [7J . Although conclusive 
evidence of exciton condensation is still lacking, two experimental groups [HI [9] have 
recently reported the observation of an upturn in the drag resistivity as the temperature 
is lowered. This upturn is interpreted as being due to strong pairing fluctuations 
that precede exciton condensation [lOj and thus serves as a precursor signal for the 
transition, similar to the enhancement of the conductivity in superconductors due to 
superconducting fluctuations above but close to the critical temperature [Hj . 

A closely related situation arises when the two layers of a two-dimensional electron- 
electron bilayer placed in a strong perpendicular magnetic fleld are close enough to 
allow the establishment of interlayer coherence [12]. In this case, the two layers in the 
system can be labelled "up" and "down" along a "z"-axis, so that the which-layer degree 
of freedom becomes a spin one-half pseudospin. Interlayer coherence in this language 
corresponds to pseudospin ferromagnetism with an easy x-y plane, since this orientation 
of the pseudospin describes a particle that is neither in the left nor in the right layer, but 
in a coherent superposition of the two. Furthermore, Coulomb drag becomes pseudospin 
drag, the mutual friction between two pseudospin states due to Coulomb scattering. 
This analogy prompted studies of spin drag, the frictional drag between electrons with 
opposite spin projection, in a single semiconductor [13j. While the realization of separate 
electric contacts to the two spin states remains an experimentally challenging problem, 
the spin drag is observed indirectly, by measuring different diffusion constants for charge 
and spin [14] . 
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Because of the presence of other relaxation mechanisms, spin-drag effects are usuaUy 
not very large in semiconductors, and are even smaller in metals. This is completely 
different in cold atomic gases where scattering between different hyperfine spins is 
the only mechanism to relax spin currents, and was considered both for fermionic 
atoms [151 ES], and for bosonic ones [TTj. In this paper we consider spin drag in a 
two-component Fermi gas, in one, two, and three dimensions. We point out that a 
particularly interesting situation occurs when spin drag is considered in a two-component 
Fermi gas that is close to a ferromagnetic instability [HI [IHl EOj ETJ [22], as can occur 
for sufficiently strong and repulsive interactions in two and three dimensions. We show 
that the spin drag is strongly enhanced as the ferromagnetic state is approached from 
the normal side [23], as expected from the analogy between electron-hole bilayers and 
pseudospin ferromagnets. In one dimension, however, where the ferromagnetic phase 
transition is absent, the effects of spin drag vanish linearly with temperature. One of 
our motivations for considering this effect is the recent observation of ferromagnetic 
correlations in a two-component Fermi gas with strong repulsive interactions [24]. The 
fact that spin-polarized domains were not directly observed adds to the theoretical 
interest [25] in this experiment. Because atoms are neutral, the relevant experimental 
quantity is the spin-drag relaxation rate, which for instance determines the damping 
rate of the spin-dipole mode in trapped cold-atom systems [16j and is thus accessible 
experimentally. Interestingly, an electronic analog of the spin-dipole mode also 
exists [26j . 

The remainder of this paper is organized as follows. We first derive an 
expression for the damping of the spin dipole mode from the Boltzmann equation. As 
mentioned before, this damping is determined by the spin-drag relaxation rate, which is 
subsequently evaluated in one, two, and three dimensions. We end with our conclusions 
and a short discussion. 



2. Spin dipole mode and spin-drag relaxation rate 

We consider a mixture of fermionic atoms of mass mmd dimensions, with two hyperfine 
states denoted by | t) and | \) . The grand-canonical Hamiltonian with external trapping 
potential V{x) and chemical potential ji is given by 



+ U jd'x i,lix)^lix)^^{x)^^{x) , (1) 

in terms of fermionic creation and annihilation operators ipl^{x) and respectively. 
At low temperatures 5- wave scattering, described by a pseudopotential V{x — x') = 
US{x — x')^ dominates, and we have therefore omitted other interaction terms from 
this Hamiltonian. We here consider only the balanced case in which there is an equal 
number of atoms in each hyperfine state. 
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Following the discussion in Ref. ^] we now derive an expression for the damping 
of the spin-dipole mode from the Boltzmann equation for the distribution function 
/cr(a3, for the atoms in spin state given by 
a/t 1 df^ hk df^ _ 
^'h^^'^^m'^- ' 
where we take the trapping potential to be harmonic V{x) = rnY^^^-^ujjx'j /2. The 
equation for is found by replacing ^ in the above. Below we give an explicit 
expression for the collision integral Fcoiif/t, /^]- 

We solve this inhomogeneous Boltzmann equation by making the ansatz 
f^{x^k^t) = NY{ek-mv^(t)/h + — ^t(^)))' with a similar expression for f^{x^k^t). 
Here, = h^k^ /2m is the single-particle dispersion and A^F(e) = [e^^^~^^ + 1]"-^ is 
the Fermi-Dirac distribution function with (3 = {k^T)~^ the inverse temperature. This 
ansatz is parameterized by the center-of-mass velocity V(j{t) and position X(j{t) of the 
atomic cloud of atoms in spin state From this, we get the equations of motion 

Nm'^ = -N'!^ + T(v,-Vi.Xi-x,); (3) 

with the function T{vi — ccj, — x-\) given by 

/r (jdj^ 
d'^x J J2^^^ 

(4) 



where 



rcoii[/,,/t] = ^^f/'/ 



(27r)'^+\,2 f d'^k2 f d'^ka f d'^k 



{2tiY J {2tiY J {2nY 
X S'^{k + ^2 - fes - k4)S{ek + Cfca - e^g - CfcJ 
X {[1 - f^{x, k, t)][l - fi{x, fea, t)]ff{x, ks, t)fi{x, k^, t) 
- f^{x,k,t)U{xM, t) [1 -h{xM, t)] [1 -h{xM, t)] } . 
We linearize the above equations using that T{v^x) V'v due to the isotropy of the 
collision integral. The linearized equations then yield a collective-mode spectrum with 
2d modes, corresponding to two types of oscillation in the rf-dimensional trap. One set 
of modes is undamped and has frequencies cjj, j G {1, . . . , rf}, and corresponds to an in- 
phase oscillation of the two clouds in the harmonic trap. The other mode corresponds to 
the out-of-phase spin dipole oscillation of the two spin states. This mode is damped as a 
result of the friction, i.e., the spin drag between the two spin states during the oscillation. 
This friction is due to collisions between particles of opposite spin and results in transfer 
of momentum between the two clouds, leading to spin drag and damping of these modes. 
These modes have the frequencies 



u;f P = + . _ ^2 . (5) 
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The imaginary part of the above frequencies gives the damping rate of the modes, and 
is given by 7 = (2rsd)~'^ = V /Nm with Tgd the spin-drag relaxation time. 

We proceed by giving an expression for Tgd for a homogeneous system with density n 
per spin state for which we, in first approximation, have taken the central density in the 
trap to make connection with the inhomogeneous case. In terms of the noninteracting 
(Lindhard) response function at nonzero temperature 

Xo{q,i^) = 2 —- ^ -, (6 

the expression for V can be worked out to yield 

rsd(T) Amnk^T J {2^Y d tt sinh^[ncj/(2A;Br)] ' 

In the next section we present results obtained by evaluating this expression. We end 
this section by noting that the above expression for the spin-drag relaxation rate is 
similar to the expression for the drag resistivity in electron-hole bilayers [6]. This also 
implies that the spin-drag relaxation rate defines a dissipative transport coefficient for 
a cold-atom system and thus represents a natural starting point for studying transport 
phenomena in these systems. 



3. Results for the spin-drag relaxation rate 

In this section we present results for the spin-drag relaxation rate l/r^d- These 
results are, among other parameters, characterized by the Fermi wave number fcp = 
[(ir(rf/2)n/47r^/^]-^/^, where r(x) is the Euler Gamma function. We also introduce the 
Fermi energy Sy = k^T^ = h^kY/'im. The results for the one-dimensional (three- 
dimensional) case are also discussed in Ref. [15] (Ref. [23j). 



3.1. One dimension 

A one-dimensional trapped gas can be experimentally realized by tightly confining two 
directions in the harmonic trap. We therefore take cji = 6^2 = (jO±_ to be much larger 
than 6^3. In the limit a <C a^, where a±_ = a/WJituJ^ and a is the three-dimensional 
5-wave scattering length, one has for the effective one-dimensional coupling constant 
that [/id = 2h^a/ma\ [28]. It is also customary to introduce the dimensionless Yang 
parameter 7 = mllm/h^n. 

In Fig. [1] we show the results that follow from Eq. ([7| by taking d = 1 and 
U = UiD- From this plot it is seen that the spin-drag relaxation rate vanishes linearly 
with temperature. It can be shown [15j from Eq. ([T]) that 

2 /cb^^ 8 2 ^ kj^T^ 



1 T^O 



rsd(T) 



-7 — + 7^7 



97r 2^F 37r \ 2^f 

in agreement with the numerical results for T/Tp <C 1 (see the inset of Fig. [T]). This 
leading order behavior is in agreement with results from bosonization [29]. The linear- 
in-T term in square brackets in Eq. (Isl) originates from contributions to the spin-drag 
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Figure 1. Spin-drag relaxation rate (in units of e-p/K) as a function of the 
temperature for a one-dimensional Fermi gas. The inset shows a zoom of the low- 
temperature region < k^T /2£y < 1, with the filled circles representing the analytical 
result in Eq. ([8|. 



relaxation rate that are controlled by momenta q of the order of 2/cf, while the quadratic- 
in-T term comes from momenta q near 0. Calculations beyond second-order perturbation 
theory have been carried out by Pustilnik et al. [30j in the context of Coulomb drag 
between quantum wires: these authors have considered only contributions to the drag 
transresistance coming from momenta q near and found p/^) oc for T ^ 0. 
Neglecting the 2/cf contributions to is fuUy justified in their case since the inter-wire 
Coulomb interaction at wave vectors of the order of 2/cf is suppressed by the exponential 
factor exp(— where I is the inter- wire distance. 

3.2. Two and three dimensions 

In two and three dimensions the spin one-half Fermi gas is predicted to undergo a 
ferromagnetic phase transition [181 UHl EQl ETJ [22], which motivated the experiments 
by Jo et al. [24]. Assuming a second-order phase transition (note, however, that the 
transition is predicted to become first order at very low temperatures [22j), the transition 
is signalled by a diverging spin susceptibility XSzsA^-,^) at zero wave vector [q = 0) 
and frequency {uj = 0). Within Stoner mean- field theory this spin susceptibility is 
calculated by summing all random-phase approximation (RPA) bubble diagrams, which 
yields xs^sAq^^) = Xo{q^^)/['^ + t^Xo(?5 ^)/2]. Hence, the critical temperature, both 
in two and three dimensions, is determined by the condition 1 + [/xo(0, 0)/2 = 0. This 
equation gives, together with the equation n = 2 J d^q A^fI^^) / (27r)^ for the total density 
that determines the chemical potential, the critical temperature Tc as a function of U. 
In two dimensions this yields 




(9) 
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Figure 2. Spin-drag relaxation rate in two dimensions as a function of temperature, 
and for various strengths of the interactions determined by ^ = Tra/a^. The vertical 
lines indicate the critical temperature for the ferromagnetic phase transition. 



with v{eY) the total density of states at the Fermi level. Note that only when 
UTy{ep)/2 > 1, the Stoner criterion, there exists a ferromagnetic phase transition. 

Experimentally, the two-dimensional situation can be achieved by tightly confining 
the system in one direction by making one (say Uz) of the three trapping frequencies 
much larger than the other two. The effective two-dimensional interaction is then 
determined by U2D = 47ra/i^/ma^, where = ^/h/mujz- 

To account for the effect of ferromagnetic fluctuations, we evaluate the effective 
scattering amplitude between atoms by summing all RPA bubble diagrams. This gives 

= U2D + ^ Xnniq.Cj) - XS.sM^Cj) , (10) 

with Xnn{Q^^) = Xo(?5^)/[l — t^Xo(?5^)/2] the RPA density response function. In what 
follows we numerically evaluate the result in Eq. ([7| with the above effective interaction, 
i.e., after making the replacement U \A^^{q^(jj)\ in Eq. 0. In Fig. |2| we show the 
results for the spin-drag relaxation rate in two dimensions, as a function of temperature 
and for various values of the dimensionless parameter ^ = vra/a^. Clearly, for sufficiently 
strong interactions, i.e., sufficiently large this rate is enhanced upon approaching the 
ferromagnetic phase transition, as discussed in the introduction. For interactions that 
do not fulfill the Stoner criterion the spin-drag relaxation rate vanishes quadratically 
with temperature, as expected for a Fermi liquid. We note that, in two dimensions, there 
is a logarithmic correction to this quadratic temperature dependence [31] although this 
is hard to discern from the numerical results in Fig. [2] and not the focus of this paper. 

The three-dimensional results are obtained by replacing U2D with the three- 
dimensional two-body T-matrix ATrati^/m [23j. The results for the spin-drag relaxation 
rate are shown in Fig. |3] and are qualitatively similar to the two-dimensional case. For 
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Figure 3. Spin-drag relaxation rate of a three-dimensional Fermi g as as a function of 
temperature T, for various values of the interaction parameter k^a. The vertical lines 
indicate the critical temperature for the ferromagnetic phase transition. 



weak interactions, such that there is no ferromagnetic phase transition, the spin-drag 
relaxation rate vanishes quadrat ically with temperature. 

4. Discussion and conclusions 

We have presented results for the spin-drag relaxation rate for a one-, two-, and three- 
dimensional two-component Fermi gas of ultracold atoms. In two and three dimensions, 
and for sufficiently strong interactions such a system may undergo a ferromagnetic 
phase transition. The spin-drag relaxation rate is strongly enhanced as this transition 
is approached from above, which could be observed experimentally as an increased 
damping of the spin dipole mode. This enhancement is determined by including 
all bubble-diagram contributions to the effective interaction between different spin 
components of the gas. This essentially treats the ferromagnetic phase transition 
within Stoner mean-field theory. In three dimensions this is most likely qualitatively 
correct, although the transition occurs at strong coupling. In two dimensions the mean- 
field results for the critical temperature are an upper bound. This is because in this 
case, and in particular in the experimentally relevant case that an external field is 
present to trap two low-field seeking hyperfine species, the phase transition is of the 
Kosterlitz-Thouless type. It is known that the Stoner mean-field theory overestimates 
the Kosterlitz-Thouless transition temperature. 

In one dimension the spin-drag relaxation rate vanishes linearly with temperature. 
In principle, we could also have included an effective interaction that included all 
bubble-like diagrams in the one-dimensional situation as well. This would result in 
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an enhancement, at some temperature Tsdw, of the spin-drag relaxation rate due to 
the divergence of the spin-density response function at zero frequency and q = 2/cf that 
signals the onset of spin-density wave antiferromagnetism. This mean-field treatment, 
however, is not accurate in one dimension, not even qualitatively. Instead, it is known 
from bosonization theory that the linear dependence at small temperatures is the correct 
one. Since this behavior is reproduced by our expression in Eq. ([7| without including 
additional contributions to the effective interactions, we do not include such fiuctuation 
corrections in one dimension. 

In future work we investigate the behaviour of the spin drag in the spontaneously 
spin-polarized phase, i.e., for temperatures T < Tc. Further studies will also investigate 
the role of critical fiuctuations close to the critical temperature, and the situation of 
negative scattering length. 
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